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ABSTRACT 

We present a first study of the star-forming compact dust condensations revealed by Herschel in the two 2° X 2° Galactic Plane fields centered at 
[I, b\ - [30°, 0°] and [i, b\ - [59°, 0°], respectively, and observed during the Science Demonstration Phase for the Herschel infrared Galactic Plane 
survey (Hi-GAL) Key-Project. Compact source catalogs extracted for the two fields in the five Hi-GAL bands (70, 160, 250, 350 and 500yL/m) were 
merged based on simple criteria of positional association and spectral energy distribution (SED) consistency into a final catalog which contains 
only coherent SEDs with counterparts in at least three adjacent Herschel bands. These final source lists contain 528 entries for the i - 30° field, 
and 444 entries for the i - 59° field. The SED coverage has been augmented with ancillary data at 24 yum and 1 . 1 mm. SED modeling for the subset 
of 318 and 101 sources (in the two fields, respectively) for which the distance is known was carried out using both a structured star/disk/envelope 
radiative transfer model and a simple isothermal grey-body. Global parameters like mass, luminosity, temperature and dust properties have been 
estimated. The L^QilM^^iy ratio spans four orders of magnitudes from values compatible with the pre-protostellar phase to embedded massive 
zero-age main sequence stars. Sources in the i - 59° field have on average lower L/M, possibly outlining an overall earlier evolutionary stage with 
respect to the sources in the i - 30° field. Many of these cores are actively forming high-mass stars, although the estimated core surface densities 
appear to be an order of magnitude below the 1 g cm~^ critical threshold for high-mass star formation. 
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^. 1. Introduction 

The chance to observe the far-infrared emission from cold dust 

condensations in star- forming regions is one of the most impor- 

. tant improve ments offered by the PACS (iPoglitsch et all 120101) 

■ and SPIRE (iGriffin et alJ , |2Q1Q|) cameras on b oard the ESA 
' Herschel Space Observatory ("Pilbr att et al.L[201Ql) . The Herschel 

Infrare d Galactic Plane Survey (Hi-GAL 'Molinari et al., 
' l2010allb) is a Herschel Key-Project that maps the inner Galaxy 
] (1^1 < 60°, IZ?! < 1°) in five photometric bands (centered at 

■ 70, 160, 250, 350, and 500 yum) and has been designed mainly 
" for observing star-forming regions and cold interstellar medium 

(ISM) structures with unprecedented spatial resolution. Thanks 
to its capabilities in terms of spectral coverage, resolution and 
sensitivity, it is expected to provide a huge progress in the study 
of the early stages of star formation, and in particular for dis- 
closing new knowledge about the formation of massive stars. 

Indeed, the mass of a young stellar object (YSO) plays a 
critical role in distinguishing between various possible regimes 
and the timescales of evolutionary paths. Whereas the sequence 
for the formation of Solar mass objects is quite well-defined, for 
high-mass protostars (M > 8 M©) it is not completely clear yet, 
due to the rarity of these sources in the Solar neighborhood (J < 
1 kpc) and their shorter evolutionary timescales (it is expected 
that such stars begin burning hydrogen while still accreting mass 
from the parental gas envelope). 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



The analysis of the infrared spectral energy distribution 
(SED) is a powerful method used to classify YSOs from 
the evolutionary point of view, and indeed the classifica- 
tion scheme (coded in Classes from to III) proposed by 
iLada & W ilking ("1984) and refined by lAndre et al.1 (Il993h is 
now well-established and commonly used for charac t erizin g the 
population of star- forming regions. ISaraceno et al.l (Il996h de- 
rived physical properties from SED fits, and plotted in a diag- 
nostic diagram the bolometric luminosity vs the envelope mass 
to easily categorize the YSOs and describ e the entire Class 0- 
II evolution. More recently Molinari et al.l (|2008) extended this 
analysis to massive YSOs, and demonstrated that it is appli- 
cable and meaningful also in this case. In this paper, we in- 
tend to apply the same methodology to the source samples ex- 
tracted from the Hi-GAL maps obtained during the Herschel 
Science Demonstration Phase (SDP), and to perform a very first 
test on its ability to describe the star formation timeline of the 
investigated regions. We describe how the SEDs were assem- 
bled and how models were fitted to data to estimate the phys- 
ical parameters of the sources. Finally, we briefly discuss the 
^boi/^env VS Mgnv diagram. 



2. SED building 

For each of the two SDP fields, a first catalog based on the 
Hi-GAL image photometry has been compiled, identifying the 
source s detected in the five diflferent bands (see Molinari et all 
l2010bl for details) based on simple positional association. 
Starting from the source extraction list at the longest wavelength. 
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Table 1. Source catalog statistics 



Stage I 
Stage 2 
Stage 3a 

Stage 3b 



^ = 30° £ = 59° 



All sources detected in at least one 
PACS/SPIREband 

Stage 1 sources excluding those with ir- 
regular SEDs 

Stage 2 sources identified in at least the 
70, 160, 250 and 350//m bands 
Stage 2 sources identified in at least 
three bands, but also with a known dis- 
tance 



2022 



528 



236 



318 



1322 
444 
141 

101 



namely 500 yum, an association with the previous band (in this 
case, 350 jim) was established if a source was found within a 
search radius corresponding to the Herschel half-power beam- 
width at the longer of the two considered wavelengths; if multi- 
ple associations were found, the closest source was kept. Other 
counterparts fulfilling the positional association criterium were 
not merged, but were added anyway as independent sources and 
considered as such in the following steps. The celestial coor- 
dinates assigned to the merged sources are those correspond- 
ing to the shorter wavelength association, which is by definition 
characterized by higher spatial resolution. After having corre- 
lated the first two wavelengths, a list composed by associated 
sources (hereafter considered as the same source) and by unasso- 
ciated single-band detections in both wavelengths was obtained. 
Entries in this list were then correlated with the next wavelength 
source list according to the same criteria, and so on, down to 
70 yum. The occurrences of multiple associations tend to be more 
frequent with decreasing wavelength (and increasing resolution), 
but for the time being we did not attempt to split up the flux at 
any given band into the contributions by difiTerent counterparts 
detected at the shorter wavelengths. 

The catalogs resulting at the end of this procedure (which we 
will call Stage 1 catalogs) contain 2022 entries for the £ = 30° 
field, and 1322 entries for the £ = 59° field. Entries in these 
catalogs are in principle SEDs exhibiting a variable number of 
counterparts for the various Herschel bands for which the po- 
sitional association was successful. An investigation of the re- 
sulting SEDs reveals a variety of situations in which the posi- 
tional association clearly did not work properly; this is apparent 
in SEDs which have missing bands and/or deep concavities due 
to a flux that is inconsistent with adjacent ones. This may hap- 
pen either because the simple positional association failed, or 
because the flux extracted for that sour ce in the detection and 
photometry stage iMolinari et al.l (1201 ObD in that particular band 
was corrupted (e.g. the source fitting did not converge because 
of source crowding or particularly difi^icult background condi- 
tions). These irregular SEDs were excluded from further analy- 
sis, as they are unreliable. The remaining sources populate the 
so-called Stage 2 catalog (see Table [T]). The subsequent selec- 
tion step is two-fold and generates on the one hand a Stage 3a 
catalog containing sources detected in at least the four 70,160, 
250 and 350 yum bands, which will be used for the color-color 
analysis. On the other hand, for the purpose of SED modeling 
(see Sect. [3]) we choose sources which have at least t hree b ands, 
but for which a distance is available (iRusseil et al.L 12010) : this 
will generate the Stage 3b catalogs. 

In Fig.[T]a [70-160] vs [250-350] two-color diagram is pre- 
sented for the Stage 3a catalog. For comparison, a grid repre- 
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Fig. 1. Two-color diagram [70-160] vs [250-350] for the Stage 
3a sources cataloged in the two investigated fields (crosses for 
£ = 30°, and diamonds for £ = 59°, respectively). Red dotted 
lines represent the loci of greybodies for various constant values 
of temperature and emissivity law exponent. In the bottom right 
corner the average error bar is reported. 



senting the isothermal greybody model loci is also overplotted. 
The analytical expression of a greybody is 



M Ko 



vo 



By(T) 



(1) 



where M is the total (gas + dust) core mass, kq is the total mass 
absorption coefiftcient evaluated at a fixed freq uency vp (in this 
case Kq = 0.005 cm^ g"^ at vq = 230 GHz, iPreibisch et al.l 
1993), d is the distance to the object, j3 is the dust emissivity 
index, and By(T) is the Planck function. 

Although there is no evidence of segregation, we can aflftrm 
that the colors are consistent with greybodies having tempera- 
tures in the range of 20 < T < 60 K. We note however that a 
noticeable fraction of sources lies in the region of the plot corre- 
sponding to negative yS values. One possibility is that the single 
temperature assumption in Eq. [T] is incorrect and that a multi- 
component fit would be warranted. This is in line with the ac- 
cepted scenario in which envelopes around YSOs show a radial 
temperature gradient. The other possibility is that we are overes- 
timating the flux at longer wavelengths (thus flattening the SED 
in the submillimeter range and mimicking a low or negative jS) 
because we are associating counterparts whose size is increas- 
ing at increasing wavelengths. Both eflfects are probably playing 
a role here. 

An alternative approach in building SEDs would have been 
to use the same size at all wavelengths, but this would have 
meant either degrading the map resolution to a common grid (the 
500yum one), or forcing the detection at 350 and 500yum using the 
sources detected shortward of 250yum and constraining the flux 
estimates to the same source size as measured at 250yL/m (and 
taking the diflTerent beam sizes into account). There are disad- 
vantages in both cases; besides, the flux extraction from the 350 
and 500yL/m maps on the locations of the 250yum detected sources 
is not at all straightforward. These alternative possibilities will 
be investigated in subsequent works. 
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3. SED fitting 

To obtain a first estimate of the physical parameters of the cata- 
log sources, the obser ved SED s were fit using both /) the grid of 
models from Robitaille et alJ |2006) and the SED-fitting tool of 
iRobitaille et alJ (20071), based on the YSO/disk/envelope model 
of I Whitney et alJ ('2003'), and //) the simple greybody model (see 
Eq.[T]). The Stage 3b source catalog (Table [TJ was used for the 
analysis: this reduced the final investigated subsample to 318 ob- 
jects for ^ = 30°, and 101 at ^ = 59°, respectively. 

In order to achieve a wider SED coverage, ancillary data 
from other Galactic plane surveys were exploited to find possi- 
ble counterparts and associate them according to the same crite- 
ria described above . In particular, on MIPSGAL maps at 24 jim 
(ICarey et al.L l2009l) a source extrac tion was performed in the 
same way as for the Hi-GAL bands (JMolinari et al.l l2010bh . re- 
turning 294 counterparts at ^ = 30°, and 89 at ^ = 59°, re- 
spectively. In addition, sources in the 1.1 mm band were re- 
trieved from the Bolocam Galactic Plane Surve y catalog (BGPS, 
iRosolowskv et al.L 120091: lAguirre et all. |2010!). A large diflTer- 
ence in number is found between the millimeter counterparts 
available in the ^ = 30° and ^ = 59° regions; the final num- 
ber of BGPS counterparts associated with Stage 3b sources in 
these two fields is 105 and only 4, respectively. This discrepancy 
arises because the £ = 59° field is relatively far from the inner 
Galaxy, compared with £ = 30°, then a smaller number of mil- 
limeter sources along the line of sight is reasonably expected. 

The model grid of IRobitaille et al.l (l2006h covers a large 
range of stellar masses (from 0.1 to 50 M©), and evolutionary 
stages (from embedded protostars to pre-main- sequence stars 
with low-mass circumstellar disks). The fitting tool uses linear 
regression to identify the best fitting models, allowing the inter- 
stellar extinction Ay and distance to be free parameters within 
user-defined ranges. All Herschel and MIPS 24 yum fluxes were 
assigned 20% uncertainties, and all models that fit with a x^ 
value satisfying x^ - xl^^t < 3 x ^data were considered good 
fits. An interstellar extinction range of to 20 was explored. 
Examples of a good and bad fits are shown in the top two pan- 
els of Fig. |2l Given the high number of free parameters in the 
models and the relatively limited number of SED points we con- 
centrated our attention on the macroscopic indications that can 
be provided by the model, namely the bolometric luminosity and 
the envelope mass. We took all values of these two parameters 
for the set of fits considered acceptable (see above) and the me- 
dian values are estimated and used as Menv and Lboi in our sub- 
sequent analysis. 

The SED-model fit provides acceptable results according to 
the above criteria for 127 source in the £ = 30° field, and 43 
sources in the ^ = 59° field. For the remaining sources no ac- 
ceptable fits were found in the entire grid of models. A visual 
inspection showed that they are mostly sources where the SED 
peaks at A >170yL/m. These situations generally correspond to 
cold envelopes with an important fraction of the gas at tempera- 
tures T <30K, not considered in the model grid of Robitaille et al. I 
(l2006l) . These SEDs were fitted with a greybody function (Eq.[T]) 
which is probably more adequate to describe early- stage (or pre- 
stellar) cores than the complex case of a protostar embedded in 
a dense envelope. The free parameters derived from weighted 
least- squares fitting (corresponding to the minimum x^ value) 
are mass, jS index and temperature. The lower panels of Fig. O 
show two examples of SED best fit from £ = 30° and ^ = 59°, 
respectively. They are chosen to display two diflferent SED ty- 
pologies: i) a source detected in all bands, with a 24 jim flux, 
which suggests an embedded warm YSO, and ii) a source whose 
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Fig. 2. Panels (a) and (b): example of good and bad fit for 
two sources from ^ = 30°, respectively, obtained with the 
IRobitaille et al.l (l2006l) fitting procedure. The black filled circles 
with error bars show the MIPS and Herschel data. The black 
solid line represents the best fit; the grey solid lines in panel (a) 
show all other models providing a good fit to the data. For the 
badly fit source only the best fit is shown. Panel (c): the same 
SED of panel (a), but this time fitted with a greybody. Panel (d): 
the SED of a source from ^ = 59°, with only SPIRE-hBGPS data, 
again fitted with a greybody. Physical parameters resulting from 
the greybody fit are also reported in panels (c) and (d), respec- 
tively. 



SED is composed only of SPIRE (250, 350 and 500 yum) and 
BGPS (1.1 mm) fluxes that can be more reliably fitted with a 
greybody model. 

As in the case shown in panel (c) of Fig.O in our sample the 
greybody fits generally fail in reproducing flu xes at wave l engths 
< 70 yum, where the approach based on .Robitaille et al.l (l2007h 
turns out to be more appropriate. 

From the fitted Stage 3b sample of sources, average temper- 
atures and dust emissivity indices were calculated. No notice- 
able diflferences emerge from the temperature distribution: the 
average temperature and standard deviation are (^30) = 23.7 K, 
o-T,o = 13.8 K, and {T59} = 22.0 K, ctj,, = 8.7 K for the two 
fields, respectively. On the other hand, the distribution ofjS val- 
ues appears to be more peaked for £ = 59° ((^59) = 1.0 ^, 
a-^3, = 0.8 K) than for £ = 30° ((/^so) = 1.6 ^, 0-^30 = 1.0 K), 
and centered on smaller values. 



4. Evolutionary timeline for massive YSOs 

The physical parameters obtained from the SED fitting can be 
used to infer the evolutionary stage of the observed sources 
by means of tools like the plot of the bolometric luminosity 
Lboi of a YSO and the total envelope mass Menv as resulting 
from the fits. It has bee n exploited by 'Saracen o et al.l (Il996|) 
and Moli nari et al.l (l2008l) to describe the evolutionary sequence 
for YSOs in the low mass and high mass regimes, respectively. 
Indeed, sources in diflferent stages are expected to occupy dif- 
ferent regions of this diagram: objects dominated by emission at 
large wavelengths (A > 250 yum) should have a Lboi /Menv ratio 
smaller than that of more evolved ones. Here we prefer to plot 
the Lboi /Menv r atio vs M^nv (Fig. O, also adapting predic tions 



obtained from iMolinari et al.l (l2008l) using the iMcKee & Tad 
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(I2OO3I) model of collapse in turbulence- supported cores. This 
model describes the standard free-fall accretion of an envelope 
onto a central core as a function of time, depicting evolution- 
ary tracks that depend on the initial value of Menv and on the 
final value of the central star mass m* when it supposedly joins 
the ZAMS. For the sake of brevity we refer the reader to the 
two papers mentioned above for a more exhaustive explana- 
tion; here we provide a brief and qualitative description. Each 
track starts from its right-bottom end, and at first proceeds al- 
most vertically. In a region of th e plane corresponding to the 
last third of the ascending tracks iMolinari et al.l (l2QQ8l) found 
sources (open grey circles in Fig.O whose SEDs were not con- 
sistent with an embedded ZAMS star, and were then proposed 
as young and accreting pre-ZAMS protostars. The end of the ac- 
cretion determines the end of the rapidly ascending tracks. In 
this region of the plot, where sources similar to YSOs associ- 
ated with UCHII regions are found, their SED is compatible with 
embedded ZAMS. Starting from here, the tracks are essentially 
determined by residual accretion and envelope dispersal due to 
molecular outflows and stellar winds. This region is populated 
with ZAMS objects whose circumstellar envelope is of lesser 
and lesser importance. 

Most of the points representing the sources of the present 
paper occupy a region Fig. [3] which corresponds to the accre- 
tion phase part of the evolutionary tracks. Crosses and diamonds 
represent sources for the ^ = 30° and ^ = 59° fields, respec- 
tively. In red are the objects which were fitted with a grey- 
bod y. In blue we plo t the s ources which are successfully fitted 
with iRobitaille et al.l (l2QQ6h models, and where the parameters 
of the central embedded forming object are not compatible with 
a ZAMS star for more than 60% of the models deemed accept- 
able. In green we plot the sources which could still be fitted as 
above, but where this time the central object has the properties 
characteristic of ZAMS stars in more than 40% of the acceptable 
models. It is reassuring that this first attempt at source classifica- 
tion confirms earlier results, in that objects presently modeled as 
embedded ZAMS are indeed located where ZAMS are predicted 
to be found. 

A quite clear distinction emerges in Fig. [3] between sources 
of ^ = 30° and i = 59°. On the one hand, an evident selection ef- 
fect is present on luminosities and masses for £ = 30° due to the 
higher typical distance (a factor 3 compared with the sources in 
the £ = 59° field), on the other hand almost all the sources from 
£ = 59° have luminosities below ~ 10^ L©, despite a wide range 
of mass values. This would suggest a global difference from the 
evolutionary point of view between the population detected in 
the two considered Hi-GAL fields, influenced by the properties 
of the star- forming regions they host. For example, evidence of 
star forn iation in a more ad vanced stage in the ^ = 30° field are 
found in lBallv et al] (l2010h . 

Figure [3] then suggests that most of the objects detected at 
present and for which the SED could be reliably determined are 
in a very early (pre-ZAMS) phase of evolution. With the present 
data it is not possible to ascertain if all detected cores are actively 
forming stars (and of which mass) or if the population fraction 
at very low L/M may be in a pre-proto stellar phase. High L/M 
ratios, however, are difl&cult to explain without an actively form- 
ing high-mass star. To be more quantitative we show in Fig.[3]the 
L/M threshold corresponding to the critical core surface density 
of 1 g cm"^ for the formation of massive stars as can be derived 
according to Krumholz & McKed (l2008l) . 

Based on this threshold, the fraction of cores actively form- 
ing massive stars would amount to 47% for £ = 30° and 13% 
for £ = 59°. However, after having calculated the surface den- 




Fig. 3. Lboi/Menv vs Menv plot for fitted SEDs in the two inves- 
tigated fields (crosses for £ = 30°, and diamonds for £ = 59°, 
respect ively). The blue and re d colors refer to fits obtained us- 
ing the IRobitaille et al.l (l2006h method (except the green points 
identifying models compatible with a central ZAMS star) and 
the greybody, respectively. The remaining grey filled circles rep- 
resent the sources studied in Molinari et al. (2008) (see their 
Fig. 9) and are plotted for comparison. Black lines represent the 
models (see text); dots mark 10^ yr inte rvals. The colored lin e 
represents the L/M threshold from .Krumholz & McKed (l2008l) 
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Fig. 4. Plot of Lboi/Menv vs core surface density in the two inves- 
tigated fields (symbols are the same as in Fig.[3l except for green 
crosses, which are not plotted here). The hashed area marks 
E = 1 g cm"^. 



sity X as the ratio between the core mass and its area in the sky at 
250 yum, we find it quite puzzling that the cores we model exhibit 
surface densities well below the critical value. Figure |4] shows 
that the surface densities calculated for most of the cores are 
lower than the critical values needed to form stars in the range 
between 10 and 200 M© (the shaded area in the figure). It can 
be argued that the adoption of different dust mass opacities kq 
could shift the distribution of points to the right, and partially 
beyond, the critical threshold. It is, however, difl&cult to explain 
why sources with L/M above the critical threshold (above the 
sloping blue line in Fig.[3l) do not occupy a specific region of the 
plot, but instead exhibit a scatter similar to the other sources in 
each region. 
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